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We examine the thermodynamic limit of fluids of hard core particles that are polydisperse in 
(N . size and shape. In addition, particles may interact magnetically. Free energy of such systems is a 

' random variable because it depends on the choice of particles. We prove that the thermodynamic 

limit exists with probability 1, and is independent of the choice of particles. Our proof applies to 
^"*^ ■ polydisperse hard-sphere fluids, colloids and ferrofluids. The existence of a thermodynamic limit 
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I . I. INTRODUCTION 

> , 

The free energy of a system of particles, as defined in statistical mechanics, depends explicitly on the the size and 
\ shape of the container holding the particles (the "system shape" ) through the dependence of the partition function 
on these quantities. Thermodynamics, in contrast, assumes a free energy density independent of the system shape 
and size. For typical systems, as the system size grows the contribution of the boundary to the statistical free energy 
becomes negligible compared with the bulk contribution. In the thermodynamic limit (system size going to infinity) 
the free energy density becomes independent of system shape and the boundary conditions. 

Ruelle |l| and Fisher M proved the existence of thermodynamic limits for a large class of fluids and solids with 
q ' interactions that fall off faster than r~ 3 at large separation. Interactions which fall off as r~ 3 or slower, complicate 
the thermodynamic limit. For neutral systems with coulombic interactions, Lieb Q proved the existence of a thermo- 
dynamic limit using the screening of the interaction at large distances. For dipolar interactions, Griffiths [Q proved 
the existence of a thermodynamic limit for dipolar lattices using magnetization reversal in a domain. The present 
i authors H extended Griffiths' proof to several non-lattice models with dipolar interactions. 

The systems discussed above contain identical particles with specified interactions. In this paper we study random 
systems, fluids consisting of hard core particles that are polydisperse in size and shape. A proof of the thermodynamic 
limit for polydisperse fluids is needed because they occur abundantly in nature and technology, and because their 
thermodynamic properties are of great interest ^ . 

For a system where the particles are chosen at random from a distribution of size and shape, the free energy is a 
random variable because it depends on the choice of particles. In general, proof of a thermodynamic limit depends 
upon the subadditivity of free energy as the system size grows. Thus we must exploit the theory of subadditive 
random variables p0[ . Furthermore, depending on the size and shape distribution of particles, a particular choice of 
particles may or may not pack into the available volume. Due to these complications, the proof requires a restriction 
on the particle size and shape distribution in the form of a "packing condition" . Thermodynamic limits of random 
systems have been investigated previously p|-^3|. Those results, however, apply only to lattice models. 

For our calculations, we fix the choice of particles in a canonical ensemble. Thus, we study "quenched" systems. 
In contrast, in an "annealed" system the relative concentration of particles of various types are controlled using 
chemical potentials The quenched and annealed free energies have been studied extensively in the context of 

KJ{ | spin glass systems p3[ , |15| |l7| . The quenched free energy is calculated by fixing the interaction between spins whereas 
the annealed free energy is calculated by averaging the partition function over the random interactions first. The 
quenched free energy is the desirable quantity since in a real system the interactions are fixed. However, it is often 
easier to calculate the annealed free energy. It is claimed that for polydisperse fluids the annealed free energy equals 
the quenched free energy Q. 
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Polydisperse fluids often phase separate into different liquid or even solid phases with different concentration of 
particle types in each phase In our terminology quenching fixes the overall concentration of various particle 

types and not theparticle concentration in each phase. This terminology should not be confused with that often used 
in the literature U of polydisperse colloidal systems to describe systems that have reached equilibrium (annealed) in 
relative particle concentration between the two phases, and systems that have not reached equilibrium (quenched). 

In this paper we prove that the free energy density F/V of any random choice of particles goes to some finite 
limiting value / with probability 1, independent of the choice of particles, as the system size goes to infinity. In Sec. || 



we carry out the proof for fluids of polydisperse hard core particles. In Sec. [II we extend the proof to fluids of hard 
core particles with magnetic interactions. These systems may be liquids such as ferrofluids Jl9[ |. Finally, in Sec. [iy| 
we summarize our results and discuss some models that are not covered by our proof. 

II. HARD CORE PARTICLES POLYDISPERSE IN SHAPE AND SIZE 

In this section we show that, for a fluid of hard core particles polydisperse in size and shape contained within a box 
of any shape, the free energy density goes to a limit with probability 1 as we take the system size to infinity. Our 
proof involves three main steps. First, we consider systems of particles contained within cubical containers and show 
that the free energy density of systems shaped as cubes goes to a limit with probability 1 as the system size goes to 
infinity. In the second step, we show that this limiting free energy density is independent of the choice of particles. In 
the the final step, we prove system shape independence of the limit by showing that the limiting free energy density 
for a system of any arbitrary shape is the same as that of a system shaped as a cube. 



In Sec. II A wc define our model that includes the particle size and shape distribution, a choice of particles from 
the distribution, the Hamiltonian, the partition function, and the free energy. We impose an artificial upper bound 
on the free energy to deal with choices of particles that do not fit inside given containers without overlapping. This 
section also defines a packing condition on the particle size and shape distribution which guarantees that with high 
probabilit y th e particles fit within containers without overlapping. 

In Sec. II B we consider systems shaped as cubes and show that the free energy density satisfies the conditions of 
a subadditive ergodic theorem ]l(J so that the free energy density of a cube-shaped system goes to a limit for almost 
every choice of particles as the cube size goes to infinity. To make use of the subadditive ergodic theorem we define a 
lattice and divide space into basic cubes defined on the lattice. A choice of particles identifies some specific particles 
associated with each basic cube. Any box defined on the lattice holds particles assigned to the basic cubes inside the 
box. The lattice, the basic cubes, and the boxes, are a device for assigning particles to containers in a random but 
well-defin ed fa shion and employing the subadditive ergodic theorems in the literature. 



In Sec. II C we show that the variance of the free energy density around its mean vanishes as the cube size goes 



to infinity, proving that the limit is in fact independent of the specific choice of particles. In Sec. [ID we prove the 
system shape independence of the free energy density by filling up an arbitrary shaped system with cubes and showing 
that the free energy density of any arbitrary shape has the same limit as that of a cube. 

A. Definitions and assumptions 

Let X be a set of particle sizes and shapes (see Fig. 1) with some probability measure /ii defined on it. Consider a 
container S, of volume Vs and a choice cjy s of Ns hard core particles contained in S, independently chosen from X, 
one at a time. Denote the set of all such choices by C/v s . The probability measure /zat s defined on Cn s is the TVs-fold 
product of the probability measure fi± defined on X. Let c denote a choice of an infinite sequence of particles and C 
denote the set of all possible choices. The probability measure [i defined on C is the infinite- fold product of 

The particles can be placed anywhere inside S with any orientation as long as the Hamiltonian Jis — T^s C 1S finite 
where 

j^hc i Q \ _ f if 110 particles overlap each other or the boundaries of S 

s ^ Ns ' 1 +oo otherwise, ^ ' 

is the hard core repulsion between the particles. Although not indicated explicitly, H s fc (cjv s ) depends on the particle 
center of mass positions {r^} and the particle orientations {fii}. Define the partition function 

i r Ns 

Zs{CNs) = m^. Ll[ dr * da ^ nsiCNs)/kBT > ® 
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where 51 = Air is the integral of dfl over all possible orientations of a particle. The prefactor iVg! is chosen to make the 
free energy extensive. This factor may be modified depending on the details of the set X to account for the entropy 
of mixing (see appendix |Aj). However, the factor N$l suffices as long as we do not "unmix" the particles 

Because the particle sizes and shapes are drawn at random, sometimes a choice of particles may not fit into the 
box S without overlapping. Some choices that fit into S may have insufficient room to move. Such choices have 
large, physically unrealistic, free energies. We wish to remove these unphysical choices from consideration. However, 
their exclusion causes mathematical inconvenience, since the probability measure for choices would no longer be the 
product of the single particle probability measure on X. On the other hand, including such choices causes the partition 
function to vanish for some cn s £ Cjv s , leading to infinite free energies. 

To handle these difficulties, we define an arbitrary threshold Zq s , for the partition function of a container S. The 
constant z a > can be interpreted as a (geometric) mean free phase space volume per particle. We define the free 
energy in the following artificial manner, 

p , s _ f -A^sfcsTlnzo if Z s (c Ns ) < Zq s ,„>. 

s[ Ns) \ -k B TlnZ s (c Ns ) otherwise. [ 1 

Provided that the particle size and shape distribution X obeys the packing condition discussed below, choices with 
Zs(cn s ) < z q S occ ur su fficiently infrequently that the limiting free energy density is independent of the arbitrary 
constant zo (see Sec. |ll C| ). 

Consider a cube T of volume Vp containing Np particles such that pr = N-p/Vr- Let Pp be the set of choices 
Cjv r £ C/v r for which Zr r (c^ r ) > z^ r , where N-p is the number of particles contained in T. Formally, 

Pr = |cAr r £ C Nr | Z r {c Nr ) > zg T X (4) 

is the set of choices that "pack" in T. Let {T} be a sequence of cubes. Our packing condition on the particle size and 
shape distribution X requires that for any zq > there exists a "critical packing density" p*(zq) > 0, such that for 
any density pr < P < P* , 

lim l x Nv {P T ) = l. (5) 

Vr — >oo 

Here pn f is the probability measure associated with the sample space Cjv r . See appendix ^ for examples of some 
particle size and shape distributions that satisfy the packing condition. 



B. Thermodynamic limit for rectangular and cubical boxes 

In this section we show that the free energy density of rectangular boxes goes to a limit with probability 1 (i.e. for 
almost every c £ C) as the system size goes to infinity. We assume that the particle distribution X obeys the packing 
condition Eq. (|B|) and that the particle density p < p* . For the proof we make use of a subadditive ergodic theorem 
(theorem 2.7 in Ref. |hJ) by Akcoglu and Krengel. 

Consider a tiling of non- negative real space by "basic" cubes with each side of length Iq. The vertices of the basic 
cubes define a lattice C such that any point in C can be written as (m, ri2, ti^Iq , where rii S Z + (nonnegative integers). 
For vectors a — (a^) and b = (bi) in £, the half-open interval [a, b) denotes the set {u\u = (ui) G C, < Ui < b{\. 
Physically, it represents a rectangular box in R^_. We call a the "base" of the box. Let I be the class of all such 
rectangular boxes. Denote the cube [w, w + 2 k d) by T™, where w is any vector in £, d is the vector (Zo, lo, lo) and k is 
any nonnegative integer. Any box / £l can be completely filled with the "basic" cubes Tq, where w runs over the 
intersection of C with /. 

Let c be a choice (infinite sequence) of particles. Break up c into successive sequences of A^o particles and put one 
such sequence into each basic cube Tq using a one-to-one map from a one-dimensional rectangular array (of sequences 
of particles) to a three dimensional rectangular array (of basic cubes) . Denote the segment of particles in Tq" by . 
For any vector u G C, define a transformation r u on C such that c' = t u (c) is another choice in C with c^ o ' = c]y u 
for every w 6 S. Since p is a product measure, t u is a measure preserving transformation on C. 

Consider a rectangular box I 6 X, of volume Vj. Let the choice of particles assigned to / be the union of choices 
c^ Q assigned to basic cubes Tq contained in / (see Fig. 2). The number of particles in / is Ni = NqVi/Vt . 
The subadditive ergodic theorem requires that for each choice c £ C the free energy Fj satisfies the following three 
conditions: 

(i) Translation invariance: For every box I £ X and every vector u £ C 
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Fi o t u = F I+U . 



(6) 



(ii) Subadditivity: For any box JeI composed of nonoverlapping boxes Ji, ...,/ n € X 

(iii) Lower bound: For some constant u>a and for all / € X 

TrFi > to A , (8) 

VI 

where Vj is the volume of box /. 
Then, the theorem states that for a sufficiently regular sequence of rectangular boxes {Ik}, of increasing size, the free 
energy density Fi k /Vi k goes to a limit /(c) almost everywhere in C as k — > oo. Note that this limit may depend on 
the choice of the particles c. A sequence of rectangular boxes {Ik} of increasing size is sufficiently regular (as defined 
by Akcoglu and Krengel flQ|0 if there exists another sequence of rectangular boxes {I' k } such that box I' k fully covers 
the box Ik, the ratio Vi k /Vji is always greater than or equal to some fixed nonzero constant, and lim^oo I' k — 

From our definition of the transformation t u , the free energy trivially satisfies (^) . For any box I £ T composed of 
nonoverlapping boxes Ii £ I, the partition functions satisfy 

Z!>Y[Z h , (9) 

i 

because each particle has more room to move inside I than an individual box f (see Ref. Q for the detailed argument). 
The free energy therefore satisfies subadditivity (]?]). The free energy also satisfies the lower bound (||) with loa — 
/csTpln(p) which is temperature times the ideal gas entropy per particle. 

Hence, the free energy satisfies all the conditions of Akcoglu and KrengePs subadditive ergodic theorem. To apply 
this to the special case of a cube, construct a sequence of cubes of increasing size {r^} such that Vr k — > oo as k — » oo 
and the distance of the bases of the cubes from the origin grows less rapidly than their size. This ensures that the 
sequence is sufficiently regular in the sense of Akcoglu and Krengel [jLM. Then, according to the theorem the free 
energy density fr k = Fr k /Vr k goes to some limit /r(c), the free energy density in a cube, for almost every ceC. 

It is interesting to note that the sequence of cubes {Tk} can be chosen so that no cubes overlap with each other 
and thus have no particles in common. This suggests that the limiting free energy density /r(c) should not depend 
on c. We demonstrate that fact in the following section. 



C. Choice independence of the limiting free energy density 

In this section we prove that the limiting free energy density /p (c) is independent of the choice of particles c. We 
show that the free energy density averaged over all choices goes to a limit and its variance goes to zero as the cube 
size goes to infinity. This implies that the limit is same for any choice with probability 1. Furthermore, the limiting 
value and its variance are independent of zq, the arbitrary constant introduced in defining free energy. 

Let (fr k ) denote the free energy density of a cube Tk averaged over all choices cpf k £ Cjv fe of Nk particles in Tk- 
From subadditivity ((?]) of the free energy it follows that 

</r fc+l > < lfr fc >. (10) 

The sequence of average free energy densities {(/r fc )} is monotonically decreasing and has a lower bound (||). The 
average free energy density therefore goes to some limit /p as k — > oo. 

To show that the variance of the free energy density of cubes vanishes as the cube size goes to infinity, we construct 
a cube Tk+i by putting together eight nonoverlapping cubes IV By using subadditivity (?]), the lower bound (@) and 
the existence of a limiting average free energy density for cubes we show (see appendix that for any e > there 
exists a sufficiently large fco such that for any k > ko, 

((/r* +1 - </r fc+1 }) 2 ) < l^fr k - (/rj) ) + e. (11) 
Iterating this procedure n times we write 
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{fr k+ j) 2 )<^((fr k -(fr k )) 2 ) + j:^e. (12) 



For a fixed value of k there exists no such that, for all n > uq, the first term on the right hand side in the inequality 
is less than e. The second term on the right hand side approaches (8/7)e from below. Hence, for any given e, there 
exist ko and no so that for any k > ko and n > Uq, the variance in the free energy density satisfies the upper bound 



2 ^ 15 



Jr k+n - (fr k+ J ) )< Y e. (13) 
For a suffi cient ly large system the variance of the free energy density becomes arbitrarily small. We previously 



proved (Sec. II B) that fr k goes to a limit /r(c) for almost every choice c of particles. A vanishing variance implies 
(by using Chebyshev's inequality |lS|]) that the probability of /r(c) differing from fr by more than some arbitrarily 
small fixed amount, also vanishes. The values of /r(c) thus converge "in probability" to /p. The existence of the limit 
fr(c) with probability 1 together with convergence in probability to fp implies that /r( c ) equals fr with probability 
1. In appendix ^ we show that fr is a convex and continuous function of the density p. 

Finally, we show that the limit fr and the variance are independent of the arbitrary constant z Q . Let Pp k be 



the set of choices that pack in a cube Tk as defined in Sec. [I A. Let Qr k be the complement of Pr k so that 



HN k (Pr k ) + MJVj,(Qrjt) = 1- Write the free energy density for a cube Tk averaged over all choices as 



fr k = / d f i Nk f Tk + / dp Nk fr k - (14) 
Jp r k JQv k 

From the definition of free energy (||) it follows that 

fr k = / dfi Nk f Fk ~ fc B Tlnz Wv fc (QrJ- (15) 

By the packing condition, PN k (Pr k ) —> 1 and ^N k (Qr k ) 0, as k — » oo. Hence fp is independent of zq. Now write 
the variance of the free energy density as 



<(/r* - fruY) = / du Nk (f Fk - fr k ) z + / d t , Nk (f Tk - f Tk Y. (16) 

Jp r k J Qr k 

From the definition of free energy (||) it follows that 

((fr k - /rj 2 ) = / dn Nk (fr k - fr k ) 2 + (fc B Tlnz + / r J 2 ^(QrJ. (17) 

Since HN k (Qr k ) ™> as k — > oo, because of the packing condition (||), the contribution of the choices cpj k S Qr k to 
the variance vanishes. Hence the variance is independent of Zq. We can then remove the free energy threshold in 
Eq. (g) by taking the limit of z a — ► 0. 



D. System shape independence of the limiting free energy density 

Now we prove that for a sufficiently regular (in the sense of Fisher Q) sequence of shapes {Sj}, the free energy 
density goes to a limit as Vsj — > oo, provided that the density of particles p = NsJVsj remains fixed. Note that the 
volumes Vs- must be adjusted so that N$- is an integer. To prove shape independence we show that the limiting free 
energy density satisfies the bounds 

lim sup f s < f r (18) 

and 



lim inf f s > f r 



(19) 



with probability 1. To show these bounds we use a technique similar to that of Fisher in his proof for identical 
particles. 

Consider a maximal filling of the shape Sj by rijk cubes T k so that each of the cubes is fully contained within 
Sj (see Fig. 3). The length Iq of the side of the basic cubes To is chosen such that there exists an integer No such 
that No/lg = p. This ensures that it is possible to achieve density p in exactly. The total number of particles in 
Sj must equal N$ j — pVsy Therefore, in general, not all rijk cubes contain equal numbers of particles. Let there be 
two types of cubes, A and B, such that type A cubes contain Njk — [Nsjnjk] particles and type B cubes contain 
Njk + 1 particles independently chosen from the particle distribution X. Here [y] denotes the greatest integer less 
than or equal to y. Note that because of our choice of Zo, ^jk must equal 2 k N plus some nonnegative integer. To 
ensure that the number of particles in Sj is Ng j , the numbers of type A and type B cubes (respectively nf k and n? k ), 
obey 



N S] = nf k N jk + nf k (N jk + 1) (20) 



and n jk = nf k + nf k . 

Let pf k = Njk/Vr k and pf k = (Nj k + 1)/Vr k be the density of particles in cubes of type A and B respectively, and 

^fk — n fk/ n jk an( i 4>fk = n fk/ n jk be the fractions of cubes of type A and B respectively. Define njkVr k /Vs j = 1 — Qk, 
so that Qk is the fraction of volume in Sj unfilled by the cubes T k . The density p of particles in Sj is related to the 
densities of particles in the cubes by 

p = (i - C]k){<j>fkpfk + fifkpfk)- ( 21 ) 

The density pf k must always be greater than or equal to p because if pf k < p, then p? k < p (since type B cubes have 
only one more particle than type A cubes), and Eq. (|2l]) cannot be satisfied for Q k > 0. In the event that Qk — we 
have <j)^ k — 1 and the equality p^ k — p holds. 

Let c be a choice (infinite sequence) of particles. Put the first Ns j particles of the sequence into Sj and denote them 
by ■ Split csj into n^ k successive sequences of Njk particles followed by nf k successive sequences of Nj k +1 particles. 
Label these sequences as c^. fc and c^-. fc+1 (depending on their length) using a map from the one-dimensional array 
(of sequences of particles) to a three-dimensional array (of cubes T™). Put particles sequences labeled by superscript 
w into cube T k . From the subadditivity of free energy it follows that 

FsM^) < £ F T ^p%^ N]k ) + Fr-(pf k ,c^ k+1 ), (22) 

weA wGB 

where the sum over w € A and w £ B includes all cubes T k of type A and type B respectively. We include the density 
as an argument for free energy and free energy density for the sake of clarity where two or more different densities 
are present in the same equation. The density is assumed to be p if not indicated explicitly. The free energy densities 
satisfy 



(23) 



fsM^) v S] < vr k (j2f^(pfk^N 3k ) + E />. : /<l<-:v..i)- 

^w£A w£B ' 

Rearranging terms we get 

fsMc Sj ) < a- ,,, u -f > /,, i/^.c; i + 4- >^ /■, i/.;;;,c! , ,2-11 



^ 3 k w£A jk weB ' 



If n? k = 0, then the sum over type B cubes on the right hand side is not present because there are no type B cubes. 

Take the limit j — > oo, holding k fixed. For sufficiently regular shapes (as defined by Fisher |Q) (jfc — > 0. Since 
pfk > P can only take discrete values, it follows from Eq. (|2l]) that there exists jo, such that for all j > jo 

pfk = P. (25) 

Therefore, <f)f k — > 1 and <j)^ k — > as j — > oo. 

The free energies /r™ (p) are identically distributed, independent random variables, because the cubes F k are 
identical in shape and the particles are independently chosen. By the Strong Law of Large Numbers |L8j 
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with probability 1 . Therefore, the limit of Eq. (|24|) becomes 



Um ±-J2fr^(fr k ) (26) 



lim sup f Sf < (fr k ) (27) 



with probability 1. Now take the limit k — > oo to get the desired upper bound (Jl8|) with probability 1. 

To show the lower bound on the limiting free energy (|l9| ) enclose the shape Sj completely inside the smallest 
possible cube Tk whose edge length is an integer multiple of 2 k l$ (see Fig 3). Fill the empty space between Sj and 
Tk with rijk cubes T k which lie completely outside Sj. Fix the density of particles in Tk at p by using two types of 
cubes T k , type A cubes that contain Nj k = [p(Vr K — ) / rijk] particles and type B cubes that contain N? k = Nj k + 1 
particles. 

The number of particles in Tk , 

Nr K = N Sj + nf k N jk + nf k (N jk + 1), (28) 
where n^ k and n^ k are the number of type A and type B cubes respectively. Dividing by Vr K on both sides gives 

V Sj . ( n jkVy k \ , , A A , ,B B \ / ooA 

P = "^r P + [—^ — ) WjkPjk + (P 3 kPjk)- (29) 

We rewrite this relation, introducing $j = VsJVr K < 1, and introducing Q k > as the ratio of volume that is 
unfilled by the cubes T k to that of Sj. Formally, (Vr K — n ]k VT k )/Vs j = 1 + Cjk- Then Eq. ( f29|) becomes 

p(l - = (1 - $,(1 + Qk))^fkPfk + ^fkPjk)- (30) 



It follows from Eq. (|3C|) that pj k > p, because if pj k < p then pq k < p and Eq. (|30|) cannot be satisfied for Q k > 0. In 
the event that Q k — then 4>^ k = 1 and pf k = p. 

Let c be a choice (infinite sequence) of particles. Denote the first Nr K particles of c by cr K . Split cr K into a 
subsequence of particles followed by n^ k successive subsequences of Njk particles and n? k successive subsequences 
of Njk + 1 particles. Denote the subsequence of N3. particles by eg. , and subsequences of Nj k and Njk + 1 particles 
by Cjy.^ and c]^. fc+1 respectively, using a map from the one-dimensional array (of sequences of particles) to a three- 
dimensional array (of cubes T k ). Put particles denoted by cs^. into Sj and particles labeled by superscripts w into 
cubes T k . From the subadditivity of free energy it follows that 

FvAP^v K ) < F s .(p,c Sj ) + Fr-(pf k ,c w Njk ) + E F Tt {pf k ,c w N . h+1 ), (31) 
Rewriting the above in terms of free energy densities and isolating fg. gives 



bl l w€A weB 



Rearranging terms further, write 

kA iB 



fsM^) > a + Qjk){h E + h E 

L we A weB J 

+ ^|/r K (p,c rK ) - E fv t {p%^ Njh ) + % E fri°(pfk,c W N jk+ i) 



(32) 



(33) 



J k weA Jk we B 

Take the limit j ^ 00 so that Q k — > 0. Since p^. > p is discrete, it follows from Eq. ([}0]) that there exists jo such 
that pj k = p for j > jo- Therefore, 4>f k ~^ 1 an d -> as j 00. For the sequence {Sj} of sufficiently regular 
shapes ||, (Vr K /Vs 3 ) > for some <5s > 0. Using the Strong Law of Large Numbers |0| we get that the inequality 

lim Mf Sj > (fr k )-Ss\fr-{fr k )\ (34) 
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is satisfied with probability 1. Now take the limit k — > oo so that (/r fc ) 
satisfied with probability 1. 



Combining the two bounds llq) and (19) we get 



I™ fs. = fr 

j-»oo 



fr- Therefore, the lower bound (fL9j) is 



(35) 



with probability 1. Therefore, the limiting free energy density / exists with probability 1 for any shape and is 
independent of the shape of the system. Its value is equal to fr . 



III. HARD CORE MAGNETIC PARTICLES POLYDISPERSE IN SIZE AND SHAPE 

In this section we extend our proof to fluids which consist of polydisperse hard core particles, as in Sec. H but m 
addition, they interact magnetically with each other. One example of such a system is a ferrofluid Jl9| , which is a 
colloidal suspension of ferromagnetic particles in a carrier liquid. The thermodynamic limit in this case is complicated 
by the 1/r 3 fall-off of the magnetic interaction where the exponent of 1/r exactly matches the dimensionality of space. 

We prove the existence of a thermodynamic limit with probability 1 for fluids of magnetic hard core particles 
using a strategy similar to that of hard core particles in Sec. |l[ We show that the free energy satisfies translation 
invariance (^), subadditivity (Q), and the lower bound (||). Once these relations are established, the rest of the proof 
of a thermodynamic limit is identical to that of the hard core particles in Sec. O and is omitted to avoid repetition. 



A. Definitions and assumptions 



Let the particles be chosen from a distribution of size, shape and magnetization, X . The magnetization is uniform 
inside each particle and has a definite relation to the orientation of the particle. The magnitude of the magnetization 
is the same for all particles. We consider both superparamagnetic particles where the magnetization Mi can reverse 



independently of the particle orientation [|20 21 fii, and non-superparamagnetic particles with magnetization fixed 
relative to particle orientation. For the non-superparamagnetic particles, we require that the particle shape have a 
two-fold rotation symmetry about any axis perpendicular to its magnetization (see Fig. 4) 

Consider a container S and a choice of particles Cjv s ■ The magnetic interaction between any two non-overlapping 
particles i,j inside S is 

U- =[ d *r [ ^(M^-Mf) _ 3(M,(r)-(r-rQ) (Mj(r') ■ (r - r')) j 
S J Vi Jv 3 I |r-r'| 3 |r-r'p J ' 

where Vi and Vj are the regions of space occupied by the magnetic material of particle i and j, and M;(r) is M; for 
r inside Vi and zero outside. Write the Hamiltonian as 

„ . , , | Ti¥ if no particles overlap each other or the boundaries of S 

n s (c Ns ) = {J 00 otherwis6j (37) 



where 



N 

Hf = £ U$ (38) 

i<j=i 

is the magnetic interaction between the particles. Define the partition function Z$ for the box in the same manner 
as in (^|). For superparamagnetic particles where the magnetization can rotate relative to the particle orientations, 
we include in f2i, in addition to the Euler angles, a discrete variable Oi = ±1 specifying that the magnetization 
is parallel (+1) or opposite (—1) to a direction fixed in the particle, and J dfli includes a sum over Oi. We may 
permit magnetization to rotate independently of the particle orientation by augmenting the integration in (0) with 
an integration over the direction of M. We require that the particle distribution satisfy the packing condition (|5|) and 
define the free energy F$ as in (||). The free energy, which depends on the choice of particles, trivially satisfies the 
translation invariance (|TJ). 
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B. Subadditivity of the free energy 

To show that the free energy satisfies subadditivity (R), we consider a container S composed of two adjacent 
nonoverlapping rectangular boxes I\ , I2 G I- Define the interaction energy between the two boxes I\ , I2 by 

7~Liij2 = — Hi! — 7~ti 2 (39) 

Let Fs(\) be the free energy of the combined system when the Hamiltonian is TLj 1 + TLi 2 plus a scaled interaction 
XHi x> i 2 . Because Fs(X) is a concave function (that is, Fg(X) < 0), it is bounded above by 

F S (X)<F S (0) + \F' S (0), (40) 

where the right side is a line tangent to the graph of Fs(X) at A = 0; here F' S (X) and Fg(X) are the first and second 
derivatives. 

As a consequence of (|40|), the free energy Fs(l) of the fully interacting system satisfies the Gibbs inequality p3j 



Fs(l)<F s (0) + [H Iul2 ] x=0 , (41) 
where Fs(0) = Fj 1 + Fj 2 is the free energy of the non-interacting subsystems, and the classical ensemble average 

[Wix,/ a ] A=o = ( 42 ) 
QN I1+ N l2NiilNi2lZiiZi J LL l J l2 f_\ 3 3 h ' h 

l — 1 J — -L 

Consider a 8 operator || which acts on any given system and has the following properties. It leaves the center of 
mass positions r; of particles unchanged, it maps the particle orientations £li on to themselves in such a way that it 
leaves the integration measure ]X dfli unchanged, and it leaves the Hamiltonian TL invariant. In addition, it reverses 
the magnetization of every particle. For superparamagnetic particles, spontaneous magnetization reversal acts as a 
9 operator. For non-superparamagnetic particles, a 180° rotation of particles about an axis of two-fold symmetry 
suffices as 9 operator. 

Applying the 9 operator to I\ but not I2 leaves both 7ii 1 and Hi 2 invariant but reverses the sign of T~Lj j 2 while 
leaving the integration measure unchanged. Hence the average in ( [42| ) vanishes and from ( f4l| ) we get 

F s <F h +F l2 . (43) 

The free energy, therefore, satisfies subadditivity (0). 



C. Lower bound on the free energy 

The lower bound (||) easily follows if the potential is stable . For stable potentials 

Us > -cjV s , (44) 

where w is some positive constant independent of the number of particles. Substituting the lower bound (Q) on Tis 
in (||) gives an upper bound on the partition function and thus a lower bound (||) on the free energy. Therefore, in 
this section we will prove the stability for our model. 

Let Ms(r) be the magnetization distribution in a box S, so that Mg(r) = M, if r is inside particle i and zero if r is 
outside any particle. Let He? (r) be the magnetic field due to the particles in S. To prove stability we make use of the 
positivity of field energy. Adding the magnetic self energy of each particle to Tig 1 gives the total energy of the whole 
system, considered as one magnetization distribution (see Ref. [p3f for identities on magnetization distributions), 

-HT^-Hf+Y^Ef^^- /d 3 rH£(r).M s (r). (45) 

i=l ^ 

Here (r) is the field, due to all particles inside the container S and 

l Ml - / d 3 rHf(r), (46) 



9 



where (r) is the field from magnetization Mj of particle i with volume Vi . 

We place a lower bound on Ti by a method similar to that of Griffiths For any magnetization distribution M(r) 
and the field H D (r) caused by it 

- - f d 3 r H D (r) • M(r) = — f d 3 r |H D (r)| 2 > 0. (47) 

2 J 87T J 

Hence 

Ns 

Hf + 2 E f S > 0- (48) 

Brown |23) rewrites the self energy in p^ ) as 

Ef f = 2ir^DfM^Ml v u (49) 

where is the demagnetizing tensor of an "equivalent ellipsoid" ; it exists for a particle of any shape [^3| , and k and 
I index the components of and M. Since is positive definite, with trace equal to 1, 

E- e lf < 2tt M 2 Vi , (50) 

so that 

N s 

Hf > -2nM 2 Y^Vi. (51) 

i=i 

If the particles are unable to fit inside S without overlap then 7is = +oo because of the hard core repulsion. Otherwise, 
if particles fit inside S without overlap, then v i — V an d 

n s = Uf > -2itM 2 V s . (52) 

The free energy of a box Jgl therefore satisfies the lower bound (||), and we have previously shown that it satisfies 
the translation invariance (||) and subadditivity ^ . The rest of the proof is identical to that of hard core particles 
polydisperse in size and shape and is omitted. 



IV. CONCLUSION 



We study the thermodynamic limit of hard core fluids consisting of particles chosen from a distribution that is 
polydisperse in size and shape. We show that a thermodynamic limit exists with probability 1 and is independent 
of the choice of particles. The existence of a thermodynamic limit implies shape independence of thermodynamic 
properties. In this section we discuss some models that are not covered by our proof, although we believe they do 
possess thermodynamic limits. 

This paper addressed fluids of particles with hard core interactions and magnetic interactions. We believe the proof 
could be extended to incorporate interactions that fall off faster than 1/r 3 . One example is a Lennard- Jones fluid 
with interaction 

A ■ Ft 

XJ. . = 1_1 _ f_i (53) 

between particles i and j, where m > n > 3 and Ay, Bij > are some constants chosen randomly for every pair of 
particles from some distribution. A "random" Stockmayer fluid is another example, where in addition to the Lennard- 
Jones interaction described above, each particle i has a randomly chosen magnetic dipole moment /i^. Another example 
is a polydisperse charged colloid in ionic solution |29|j3C|] with an interaction such as 



u = mi^. (54 ) 

er 



where n > is the Debye screening length. 
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For such systems the free energies for a container S composed of two adjacent nonoverlapping rectangular boxes 
I\ , I2 el separated by a distance R > Rq > satisfy 

F s < F Tl + F l2 + Ai 2 . (55) 

Here 

N h N l2 5 12 

with some e > and <5 12 depending upon the probability distributions of the pair potentials. Our proof of thermo- 
dynamic limit uses a subadditive ergodic theorem which requires strict subadditivity (Q) with A12 = 0. However, 
we suppose that provided that the probability distributions are sufficiently restricted, a subadditive ergodic theorem 
similar to that of Akcoglu and Krengel [[l0| but with the weaker subadditivity condition ( |55| ) will hold and a proof of 
a thermodynamic limit will follow. 

Granular magnetic solids provide other examples of systems outside the realm of our present considerations. Con- 
sider a solid system with particles frozen in space but magnetic moments able to rotate with respect to the body of 
the particles [^7| (or a solid containing rotationally free magnetic particles inside cavities Q) One way to prepare 
such systems is to freeze the carrier fluid in a colloidal suspension of magnetic hard core particles so that the positions 
of particles are representative of an equilibrium configuration of the fluid. Our proof does not directly apply to such 
systems because the infinite system cannot be constructed from finite cubes of frozen particles. A system constructed 
from finite cubes will have no particles overlapping the boundaries of these cubes and thus the spatial arrangement of 
particles differs from an equilibrium configuration of the fluid. If the particles are "unfrozen" and frozen back again 
after their positions have relaxed, then subadditivity (^) no longer holds strictly, though it will hold on average. 
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APPENDIX A: EXAMPLES OF ALLOWED PARTICLE DISTRIBUTIONS 

Here we present some examples of particle size and shape distributions which obey the packing condition. Our 
discussion hinges on the ability to pack a monodisperse collection of hard spheres with a hexagonal close packing 
fraction of <p hcp — 7^7^- Our presentation is informal and primarily intended to motivate the form of our packing 

condition, rather than to achieve rigorous proofs. 

Binary mixture of hard core spheres: Consider spherical hard cores of two types with radii r a and r@ occurring 
with probabilities p a and pp, respectively. A given choice cjy G Cn of N particles yields N a spheres of type a and 
N/3 = N — N a of type f3. Define the partition function in a volume V as 

Zv (c„) = -±- [Udne-^y^. (Al) 



NJNplJv 

The factor of N a \Np\ is chosen instead of N\ because it yields the usual entropy of mixing [24[] . Given a choice cjv, a 
volume U(cjv) exists that is sufficiently large that the partition function exceeds a threshold 

Z v (c N ) > z?. (A2) 

We construct the volume V(c/v) satisfying the inequality ( |A2| ) by assigning particles to spherical shells of radius 
r' a = r a + ro and r'p = rp + ro where ^-Tq = Zq. Shells of type a and ft surround volumes v' a = and v'p = 

If we choose the volume 

V(c N ) = (N a v' a + N v' )/<p hcp (A3) 



then the inequality (A2) holds 



To verify this claim, divide the volume into subregions of volumes N a v' a /(f> hcp and N pv'p / <p hcv respectively. Closely 



pack the a- and /3-type shells into their respective subregions and evaluate the configurational integral in Eq. (Al) in 
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two stages. First, confine each particle within its shell and integrate, yielding a contribution of Zq . Second, permute 
identical particles within identical shells, yielding a factor of N a \Np\. Since the combinatorial factor cancels against 
the normalization of the partition function in Eq. (Al), the partition function Zy(c]y) is at least as large as Zq 



Now we show that, with probability 1, V(cn) does not greatly exceed Nv'(zo)/(p hcp , where we identify the mean 
volume per shell v'(zq) = p a v' a + ppv'p. Note that v'(zq) approaches the mean particle volume v as zq — > 0. Take a 
cube r of volume 

where e is an arbitrarily small positive number. Divide the set Cn of choices cat into a "packing subset" 

P T = {c N e C N | V(c N ) < V r } (A5) 

and a complementary "nonpacking subset" 

Qr = {cat e C N | V(c N ) > V r }. (A6) 

Since Pr and Qr are complementary subsets, their measures obey /i_/v(-Pr) + fi>N(Qr) — 1- By Chebyshev's inequal- 
ity []l8f on the sum of independent random variables we bound the measure of the nonpacking subset by 

^ (0r) ^ (A7) 

For sufficiently large N the probability that a choice does not pack into Vr vanishes. Thus V(cn) < Vr with probability 
1 at density 

(hhcp 

p = N V r = - \ - . (A8) 
v'{z ){l + e) 

Therefore, binary mixture of hard spheres obeys the packing condition (|b|) with p*(z ) = (f) hcp /v'(z a ). Note that 
p*(zo)v approaches (f) hcp as zo — * 0. 

Finitely many types of hard spheres: Consider some finite number, t, of types of hard spheres, each type 
with a different finite radius and a certain probability of occurrence. The argument presented above for a mixture of 



two types of spheres generalizes easily to the case of t > 2. The combinatorial factor in the partition function (Al) 
generalizes to n^iVc,! and the mean volume per shell generalizes to v'(zq) = ^2 a p a v' a - An argument similar to that 
for a binary mixture of hard spheres shows that this distribution satisfies the packing condition (g) with a critical 
packing density p*(zq) = <f) hcp / v' (z$) . A packing fraction of <fi hcp is achievable in the limit zq — ► 0, but for finitely 
many types of hard spheres a high er p acking fraction should be achievable because smaller particles can occupy the 
interstitial sites of larger particles |gf| . 

The distributions above satisfy the packing condition using the prefactor H a N a l to define the partition function. 
However, the prefactor of (%2 a N a )\ = Nl may also be used to define the partition function as in Eq. (0). The 
logarithm of the ratio of the two prefactors is the entropy of mixing between the particles of various types. As long 
as the particles are chosen randomly from the distribution and not "unmixed", the factor of N\ suffices to give an 
extensive entropy P,p4|- To show that the distribution satisfies the packing condition using prefactor TV!, enclose 
particles in shells so that they have free volume of at least z' inside the shell. The partition function is bounded from 
below by 

^>^>[C (A9) 



TV! u " V t , 

Set Zq = tzQ to show that the distribution satisfies the packing condition (||) as long as t is finite. Now the critical 
packing density is p*(zo) = <t> hcp /v'(tzo). As before, p*(zo)v approaches (f> hcp as zq — > 0. 

Continuous distribution of hard spheres with an upper size cut-off: Consider a continuous distribution 
of hard spheres with volumes < v < v ma x ■ This distribution has essentially an infinite number of types of particles 
and so our discussion above does not directly apply to it, however, it still satisfies the packing condition (g). To show 
this, define the partition function with prefactor Nl to make the free energy extensive. Break the distribution into a 
finite number, t, of bins with boundaries < v\ < i>2... < vt — v max . Let the i th bin Bi contain all the particles with 
Vi-i < v < Vi. Enclose every particle in bin Bi within a shell of volume 
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3 I V 47T 7 V 47T 



(A10) 



so that the largest particle in Bi has a free volume of at least z' Q = tzo, then apply the same argument as for finitely 
many types of spheres. The resulting p*(z ) = <f> hcp /v'(tz ) results in a suboptimal achievable packing fraction 
4> = p*(zq)v = (f> hcp v/v'(tzo) < 4> hcp due to the bin width and the choice of prefactor Nl. By using a sufficiently large 
but finite number of bins and taking zq — > 0, a packing fraction arbitrarily close to <fr p can be achieved. As in the 
case of discrete distributions, for a continuous distribution of hard spheres one expects a higher packing fraction than 
that of a monodisperse system should be achievable |^5| . 

Continuous distribution of hard spheres with no particle size cut-off: Now wc consider a continuous 
distribution of hard core spheres such that the particle size has no upper limit. To ensure that the large sized particles 
occur sufficiently infrequently, we assume that there exists a particle volume v* such that the probability of a particle 
with volume greater than vq for any vq > v* falls off as 

> «b) < -i^y. ( A11 ) 

where 7 and 5 are positive constants. Note that this condition guarantees the existence of a mean particle volume 
but not the existence of higher moments. 

To show that this distribution satisfies the packing condition (||) we break the distribution into a finite number, 
t + 1, of bins with boundaries < v\ < v^-- < Vt-i < v t < 00. We presume that Vt-i > v* . The next-to-last bin Bt 
includes all the particles with Vt-i < v < v t and the last bin Bt+i contains all the particles with vt < v < 00. We 
now send Vt to infinity in such a manner that Bt+i is empty with high probability but the maximum particle size in 
Bt grows in a controlled fashion. To show that the distribution satisfies the packing condition (||) we first show that 
the set of choices with no particle in Bt+i is a set of measure 1. Then we show that in this set there is a subset of 
measure 1 of choices which have a partition function of at least Zq . 

The probability that a collection of N particles randomly chosen from the distribution contains at least one particle 
in the last bin Bt+\ is bounded by 

Nj 

Pn{cn I at least one particle in B t +i) < 1+s ■ (A12) 

v t 

We wish to avoid all choices with at least one particle in Bt+i- To ensure that the probability of such choices goes to 
zero for large N, Vt must grow faster than N 1 ^ 1+s \ At the same time v t must grow slower than N to ensure that 
the largest particles in Bt are much smaller than the system size. Hence let vt grow proportional to N 1 ^ 1+ ^'> , so that 

7' 

Pn(cn I at least one particle in B t +i) < 5 ► 0, (A13) 

where 7' is some constant. For large N, the probability approaches 1 that no particle is in bin Bt+i- 

Now consider the particles in bin B t ■ The next-to-last bin B t must be treated specially because the largest particle 
size in Bt diverges as the system size to goes infinity. Enclose each particle j ( with v olume Vj) in Bt inside its own 
shell with volume v'j so that it has a free volume of z' Q — tzo. Using condition ( All ) the total volume of the shells 
corresponding to bin B t (averaged over choices c^r) obeys 



' t )<-^^- (AH) 



V 3 

VjEBt 

Since the volume of the system is proportional to N, the volume fraction occupied by shells corresponding to bin B t 
is proportional to v^\. This volume fraction can be made arbitrarily small by choosing a sufficiently large Vt-\. We 
assume that the N t particles in bin Bt (enclosed in their shells) can be packed inside a small region of the system in 
Ntl different ways with high probability (tending to 1 as N goes to infinity). 

The particles in the remaining bins B\, Bt-i have a partition function of at least z ~ * in the remaining volume 
of the system with high probability (tending to 1) as we showed previously for a continuous distribution of spheres 
with an upper size cut-off. Therefore, these particles together with the particles in Bt satisfy the packing condition (j|). 

Non-spherical particles: Finally, consider a distribution of non-spherical particles. To show that the distribution 
satisfies the packing condition (^), enclose each particle inside a sphe rical shell with a diameter equal to the largest 



dimension of the particle. If the volume of the shells satisfy condition ( All ) then the rest of the argument is the same 
as that of a continuous distribution of hard spheres. Although this argument shows the existence of a finite particle 
number density p*(zo), it does not guarantee the existence of a finite packing fraction. In fact, it is possible to have 
distributions with oddly shaped particles (fractals for example) so that a finite p*(zo) exists but the particle packing 
fraction is zero. 
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APPENDIX B: VARIANCE OF FREE ENERGY DENSITY FOR CUBES 



To show the variance relation (|ll| ) consider two adjacent nonoverlapping cubes T l k and T™ such that 1^ — T l k U r™ 
is a rectangular box. Using the subadditive relation (^) we get 

flu < \ifrl + frf). (Bl) 



Adding ui A (see lower bound (0J) to both sides of (Bl) and squaring gives 



ff k + 2u A f Ik +u\< i(/ r , + / rr ) 2 + a£ + io A (f rl + fr T ). (B2) 



Averaging both sides over all choices of particles we write 

(fl) + 2Lo A (f Ik )+Lo\ < l(fl) + \(f Tk ) 2 + 2u A (f Th ) (B3) 
Subtracting lo 2 a + (ff) + 2tu A (fj k ) and rearranging terms on the right-hand side gives 

((fi„ ~ (//J) 2 ) < \i(fr k </r,.)) 2 ) + ifr k ) 2 (fi k ) 2 + 2w A ((fr k ) - (f Ik )). (B4) 

Using the bound 

</r fc > + (fh) < -2fc B Tlnz = 2lu b , (B5) 

we get 

((fi k ~ (fi k )) 2 ) < \((fr k - (/rj) 2 ) + 2^ A +uj B )((f Tk ) (f Ik )). (B6) 
Both (fi k ) and (Av ) go to the same limit /r since (fr k+1 ) < (// fc ) < (fr k )- Therefore, we can bound the second term 



on the right in (B6) below any e' > for sufficiently large k. It follows that 



{{fi k -(fi k )) 2 )<\{{fv k -{fv k )) 2 ) + e'. (B7) 

By stacking two rectangular boxes Ik together to form a square slab, then putting together two adjacent slabs, we 
construct a cube The variance of the cube T^+i is related to that of r\. by 

{(fr k+1 - </r fc+1 )) 2 ) < \{{fr k - {fr k )) 2 ) + e, (B8) 

where e = (1 + \ + |)e' = \e' . 

APPENDIX C: CONVEXITY AND CONTINUITY OF THE FREE ENERGY DENSITY FOR CUBES 

Consider a cube with volume Vr k that contains AT fe particles. Since is an integer, the free energy density 
fr k {p,C]y k ) is defined only for discrete values of p. Since fr k (p,CN k ) is a random variable, a linear interpolation to 
include fractional number of particles is not possible. However, the average free energy density (fr k (p)) is not random 
and we define (/(/£>)) f° r a U values of p by linear interpolation between values of p that correspond to integer values 
ofiV fe . 

Now consider a cube T^+i consisting of 4 cubes Tk each containing NL particles and 4 cubes each containing 
N k particles. From subadditivity (^) of free energy, and averaging over all possible choices of particles in the cubes 
Tk, it follows that 

^ ^((/r fe (p')) + (/r fc (/))), (CI) 



where p' = N' k /Vr k , p" — N k /Vr h and p — (p' + p")/2. By linear interpolation Eq. (CI) holds for any p' and p" 
Taking the limit k — > oo gives 
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/r(p))<^((/r(p')) + (/r(p")))- 



(C2) 



Therefore, (fr(p)) = fv{p) is a convex function of p. Since fr(p) is a bounded function of p, it is also continuous 
function of p [Em. 
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FIG. 1. An example of a size and distribution of particles. The hard core particles A, B, C are chosen with fixed probabilities. 
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FIG. 2. One particular choice of particles from the distribution in Fig 1. into basic cubes, with iVo = 2 particles in each 
basic cube. The box / =[(lo, lo), (3io , 2£o)) m this example is assigned particles C, B, B, and A. For easy visualization we show 
the figure in two-dimensions instead of the actual case of three dimensions. 
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FIG. 3. (a) A cross-section of an arbitrary shape Sj. The space inside Sj, and the space between Sj and T K , is maximally 
filled by cubes Tk ■ Dashed lines indicate the cubes that intersect the boundary of Sj . The volume of cubes Tk that intersect 
Sj becomes negligible in comparison to the volume of Sj as the size of Sj goes to infinity and the size of I\ grows less rapidly 
than Sj. (b) A cube T k in d dimensions consists of 2 kd basic cubes Fq. In this case k — 1 and d — 2. 




FIG. 4. (a) Particle shapes with a two-fold symmetry perpendicular to their magnetization. The arrow indicates the direction 
of magnetization, (b) Particle shapes lacking two-fold symmetry perpendicular to their magnetization. 
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